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Abstract 

The notion that microparsec-scale black holes can be used to probe gigaparsec-scale physics may seem counterintuitive, 
at first. Yet, the gravitational observatory LISA will detect cosmologically-distant coalescing pairs of massive black 
holes, accurately measure their luminosity distance and help identify an electromagnetic counterpart or a host galaxy. 
A wide variety of new black hole studies and a gravitational version of Hubble's diagram become possible if host 
galaxies are successfully identified. Furthermore, if dark energy is a manifestation of large-scale modified gravity, 
deviations from general relativistic expectations could become apparent in a gravitational signal propagated over 
cosmological scales, especially when compared to the electromagnetic signal from a same source. Finally, since inspirals 
of white dwarfs into massive black holes at cosmological distances may permit pre-merger localizations, we suggest 
that careful monitoring of these events and any associated electromagnetic counterpart could lead to high-precision 
cosmological measurements with LISA. 
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1. Introduction 

Essentially all astronomical measurements are 
performed via electromagnetic waves. The availabil- 
ity of accurate gravitational wave measurements 
within the next decade or so will thus be a sig- 
nificant development for astronomy. In particular, 
since the propagation of photons and gravitons 
could differ at a fundamental level, gravitational 
waves emitted by cosmologically-distant "space- 
time sirens," such as coalescing pairs of massive 



1 K. Menou gratefully acknowledges more than a decade of 
exchanges with Jean-Pierre Lasota on the physics of black 
holes (and definitely white dwarfs). 



black holes, could be used as valuable new probes 
of physics on cosmological scales. 

Black holes with masses > 1O 6 M are present 
at the center of nume r ous n e arby galaxi e s (e.g . 



Kormendv fc Richstond . Il995c iMagorrianl 119981 ) 



As such galaxies collide over cosmic times, their 
central black holes coalesce, releasing > 10 58 ergs 
of binding energy in the form of gravitational waves 
(hereafter GWs). To measure the GWs emitted 
by these cosmologically-distant space-time sirens, 
ESA and NASA will build the Laser Interferometer 
Space Antenna, LISAT^I. 

GWs emitted by black hole binaries have the unfa- 
miliar property of providing a direct measure of the 
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Fig. 1. Orbit of the LISA 3-base space interferometer, trail- 
ing the Earth around the Sun. A typical coalescing black 
hole pair observed for a year is well-modulated by the inter- 
ferometer and thus accurately localized on the sky. [Credit: 
http : //lisa . nasa . g ov/] 



luminosity distance, Dl, to the black holes, without 
extrinsic calibration. Ow ing to the hig hly coherent 
nature of GW emission (jSchut j . Il986l ), the ampli- 
tude (or strain), h +x , frequency, /, and frequency 
derivative, /, of the leading order (quadrupolar) GW 
inspiral signal scale as 



h +x (t) cc 



[(l + z)Af c ] 5/3 / 2/3 



/(t)cx[(l + z)Af c ] 5/3 / 11 /3 j 



(1) 
(2) 



where + x represents the two transverse GW polar- 
izations, M c = (mim2) 3 '' 5 /( m i+ ffl 2) 1 ' 5 is the black 
hole pair "chirp" mass and z its redshift. Provided 
the GW source can be reasonably well localized on 
the sky, an extended observation of the chirping sig- 
nal leads to precise measurements of h +x , /, / and 
thus Dl, independently. As illustrated in Fig. [TJ 
LISA's orbital configuration allows for a "triangu- 
lation" of GW sources on th e sky, to with i n a solid 
angle £Q ~ 1 deg 2 typically (jCutlerl , Il998t IVecchiol 
20041 ). This permits very accurate measurements, 
e.g. distanc es with errors 5Dl/Dl < 1% at z < 
2 typ ically (I Cutler! 1 19981: iHughesl . 120021: IVecchiol 
2004 : lLang fe Hushed . l2006lT Masses are indepen- 
dently dete rmined to ver y high accuracy (typically 
< 1%; e.g.. lHughedl2002h 



2. Post- and Pre-Merger Localizations 

In principle, the same sky localization that helps 
determine the distance to a source accurately can 
be used to find the host galaxy of a pair of merging 
black holes seen by LISA. The secure identification 
of the host galaxy would enable a wide variety of 
new galactic black hole studies (see £)3.ip . 

Initially, the prospects for finding the host galaxy 
of a pair of merging black holes were considered 
to be poor, simply because of the large number 
of galactic candidates loca ted in t h e ~ 1 deg 2 
LISA sky error-box (e.g., ICutlerl . 1998; IVecchiol 
20041) Recently, however, this possibility has been 
reconsidered, with more optimist ic conclusi ons 
dHolz fe Hughesl 120051 : iKocsis et al!l200ril2007l) . 

Given a cosmology, it is possible to translate the 
accurate luminosity distance measurement to the 
GW source into a narrow redshif t slice in which 
the host galaxy must be located ( Holz fe Hughesl 
20051 : iKocsis et al! l200d) ■ Various contributions to 
the redshift errors that arise in performing this con- 
version are shown in Fig. [2 for a representative 
equal-m ass binary, as a fun ction of the GW source 
redshift (|Kocsis et al!l2006f ). At redshifts z > 0.25, 
where most black hole binary sources are expected 
to be found, weak lensing errors due to line-of-sight 
inhomogeneities (on top of the smooth average cos- 
mology) are the main limitation to an accurate de- 
termination of the redshift slice in which the host 
g alaxy ought to be loc ated. 



Kocsis et al 



(|2006l ) have studied in detail the 
possibility that the three-dimensional information 
available (sky localization + redshift slice) could be 
used to single out a quasar, or any other unusually 
rare object (such as a star-bust galaxy), in the LISA 
error box, after coalescence. Finding such a statisti- 
cally rare object post-merger would make it a good 
host galaxy candidate for the newly-coalesced pair 
of black holes. 

However, it maybe much more advantageous to 
use a pre-merger strategy to identify the host galaxy 
of a pair of coalescing black holes seen by LISA. In- 
deed, one can use near real-time GW information 
on the sky localization, in combination with the ac- 
curate timing of the inspiral event, to predetermine 
well in advance where on the sky the merger is lo- 
cated. A unique host galaxy identification could then 
proceed through coordinated observations with tra- 
ditional telescopes, by monitoring in real time the 
sky area for unusual electromagnetic emission, as 
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Fig. 2. Contributions to the error on the inferred redshift of 
an electromagnetic counterpart to a LISA coalescence event, 
as a function of redshift z, for mi = m.2 = 1O 6 M0. The in- 
trinsic LISA error on the luminosity distance, d^, is shown 
for two representative cases (a & b, solid lines). Errors due 
to the peculiar velocity of the source (for v = 500 km s _1 ; 
short— dashed line), uncertainties on the background cosmol- 
ogy (long— dashed line), and errors due to weak lcnsing mag- 
nification (dash-dotted line) are also shown (see lKocsis et all 
120061 for details). 



the coalescence proceeds. 

A variety of mechanisms exist through which 
disturbed gas in the vicinity of black hole pairs 
will power electro magnetic emission duri n g and 



after coalescence (lArmitage fc Nataraian , 
Milosavlievic fe Phinnev . 20051: Dotti et al 



2002 



2006 



Bode fe Phinnevl . l2007tlMacFadven fe Milosavljevid 
2008). For example, at the time of coalescence, 
> 10 53 ergs of kinetic energy are delivered to 
the recoiling black hole remnant and its environ- 
ment, for ty p ical r e coil ve l ocities > 100 km/ s (e.g 
Baker et all l2006|. [iool ICampanelli et all 1200 



20071: 
annd. 



Herrmann et al. . 20071 : Schnittman fc Buonann 
2007) . This may lead to detectable signatures 
(jLippai et al. , 2008; Mi losavlievic fc Phinnevl . l2005l) 



and permit the coincident identification of a unique 
host galaxy. The detailed nature of such electro- 
magnetic counterparts remains largely unknown, 
however. 

To a large extent, LISA's ability to localize a 
long-lived source on the sky is related to the GW 
signal being modulated as a result of the detec- 
tor's revolution and change of orientation when 
the constellation orbits around the Sun (Fig. [T]). 
Even though most of the GW SNR accumulates 
during the final stages of inspiral/coalescence for 
typical GW sources, reasonably good information 
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Fig. 3. Evolution with pre-merger look-back time, tf-tiscOi 
of LISA source localization errors, for M = 2 X 10 6 Mq and 
z = 1. The top panel shows luminosity distance errors and 
the bottom panel shows sky position angular errors (equiv- 
alent diameter, 2\/afe, of the error ellipsoid). Best, typical, 
and worst cases for random orientation events represent the 
10%, 50%, and 90% levels of cumulative error distributions, 
respectively. Errors for worst case events effectively stop im- 
proving at a finite time before merger, even though the sig- 
nal-to-noise ratio accumulates quickly at late times. Errors 
for best case events (especially the minor axis) follow the 
sign al-to-noise ratio un til the final few hours before merger 
(see lKocsis et alj|2007l for details) 

on sky localizations must be available well before 
final coalescence since this information accumulates 
slowly, over the long signal modulation (orbital) 
timescale. Because of significant cross-correlations 
between sky localization and distance errors, it 
turns out that this argument is also largely valid for 
luminosity distance erro rs (jHolz fc HughesL l2005t 
iKocsiset alll2006Ll2007l) . 

Figure [3] shows the pre-merger time evolution of 
luminosity distance and angular sky localization er- 
rors for a representative black hole pair at z = 1. 
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Errors improve quickly at early times but their evo- 
lution slows down considerably at late times. Ac- 
cording to both panels, even accounting for random 
orientations of various source and detector angles 
(shown as best, typical and worst cases), signifi- 
cant information is a vailable days to wee ks prior to 
the final coalescence ( Kocsis et al. , 20071) . Including 
black hole spins in the analysis has been shown to 
result in significant improvements on the errors dur- 
i ng the last few days t o hours prior to coalescence 
(|Lang fc HughesLl2008h . 

With the expected availability, by the time LISA 
is operational, of sensitive large field-of-view (FOV) 
astronomical instruments for weak lensing and su- 
pernova studies, it becomes interesting to estimate 
the amount of time prior to merger during which 
the LISA sky localization falls within the FOV of 
such an instrument. When this happens, continuous 
monitoring of the designated sky area, unt i l fina l 
coalescence, becomes possible. iKocsis et al. ( 2007 ) 
have performed a detailed analysis of this possibil- 
ity, using LSST and its 10 deg 2 FOV as a reference. 
Figure 3] shows results for representative equal-mass 
binaries, as a function of their total mass and red- 
shift. The various contours show that prospects for 
electromagnetic monitoring days to weeks before the 
coalescence are good for sources at redshifts z < 2 
3. Monitoring for the bes t GW sources out t o z ~ 
5-7 may even be possible (jKocsis et all 120071 ). 




3. New Science with Electromagnetic 
Counterparts 



Fig. 4. Contours of advance warning times in the total mass 
(M) and redshift (2) plane for equal-mass black hole binaries 
(mi/7712 = 1). The contours trace the look— back times at 
which the equivalent radius (2\/a6) of the LISA localization 
error ellipsoid first reaches an LSST-equivalent field-of-view 
(3.57°). These contours correspond to the 50% level of cu- 
mulative distributions for random orientation events. The 
contours are logarithmically spac ed in days and 10 days is 
highlighted with a thick line (see IKocsis et al.|[2008l for de- 
tails). 

us about the nature (dry/wet) and the timing of 
the galactic merger that resulted in the black hole 
binary coalescence. Finally, measuring velocity dis- 
persions, <r, for several host galaxies, together with 
the black hole masses known from the GW signal, 
would allow us to accurately map the evolution of 
the Mbh-c relation with cosmic time, at least for 
such transitional objects as the hosts of coalescing 
black hole pairs. 



3.1. Galactic Black Hole Astrophysics 



3.2. Gravitational Hubble Diagram 



A large variety of new galactic black hole astro- 
physics would be enabled by successful identifica- 
tions of the host galaxies of coalescing black hole 
pairs. We mention only a fe w possibilities here and 
refer the inte r ested reader to lKocsis et al.1 (|2006l ) and 
Kocsis et al.l (|2008l ) for additional discussions. 

From the black hole masses, spins and binary 
orientation, all accurately constrained by the GW 
signal, one would be able to study the physics of the 
post- merger accretion flow onto the r e mnant black 
hole (jMilosavljevic fc Phinnevl . 120051 : iDotti et af 



20061 ) with unprecedented accuracy. This would 
include precise constraints on the Eddington ratio 
of the accreting source, its emission and absorp- 
tion geometries and possibly its jet phenomenology. 
Similarly, studies of the galactic host might tell 



Another consequence of successfully identifying 
the host galaxies of coalescing black hole pairs is the 
possibility to draw a gravitational Hubble diagram, 
i.e. one that relates the gravitational luminosity dis- 
tances, Dl, of these GW sources to the electromag- 
netic redshifts, z, of their host galaxies. 

One of the main interests of a gravitational Hub- 
ble diagram arises from its immunity to common 
systematics affecting electromagnetic measure- 
ments. Indeed, a gravitational Hubble diagram, 
which is based on gravitational distance measure- 
ments with self-calibrated sources, is not susceptible 
to any significant bias from absorption, scattering 
or reddening of GWs. 

In practice, however, the value of such a dia- 
gram is limited by line-of-sight matter inhomo- 
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geneities, which generate weak lensing uncertain- 
ties on the gravitational Dp measurement (see 
Fig. [21 iHolz fc Hughes!. 120051: iKocsis et al l. 12006 



Dalai et al. . 2006t iGunnarsson et al. . 2006 : iLinder . 
20081 ). While the lensing efFect can in principle be 
averaged out over many random lines-of-sights, 
it may not be possible to do so for coalescing 
pairs of massive black holes if LISA merger event 
rates are modest (e. g., a few tens per y ear at 



z < 5: iMenou et al I. l200lUWyithe fc Loebl . 12003 
Sesana et all 120041 : iMicic et all 120071 ). Weak lens 
ing errors on individual measurements amount to 
distance uncertainties ranging from 5Dl/Dl — 1% 
at z = 0.5 to SD tJDt, ~ 10% at z = 5 (e.g. 
Kocsis et al. . 2006), which makes a gravitational 
Hubble diagram imprecise even at moderate red- 
shifts. The extent to which LISA events can be used 
to draw a meaningful Hubble diagram will thus de- 
pend strongly on the actual distribution of massive 
black hole merger events with redshifts and the cor- 
responding efficiency of host galaxy identifications. 
As we describe in £|4j however, white dwarf spiral- 
ing into massive black holes may offer a practical 
avenue for precision cosmology with LISA. 



3.3. Diagnostics of Modified Gravity 

The possibility that the accelerated expansion of 
the Universe results from a failure of general rel- 
ativity has fueled much theoretical work on large 
scale modifications of gravity over the past few years. 
Since building a satisfactory theory of modified rela- 
tivistic gravity is a formidable task, any insight that 
can be gained from direct observational constraints 
on the linearized GW regime cannot be overlooked. 
LISA, with its ability to measure the GW signal from 
cosmologically-distant sources, may thus be one of 
our be st probes of modified grav ity on cosmological 
scales (jPeffavet fc Menoul . [20071) . 

One may expect gravity modifications to contain 
a new length scale, R c , beyond which gravity de- 
viates from general relativity. In order to explain 
the observed accelerated expansion of the Universe, 
this scale is expected to be of the order of the cur- 



rent Hubble radius, . An existence proof of 
modifications of this type is g iven by DGP gravity 
(jDvali et all l2000t iLud . 120061 ) , a braneworld model 
wi th an infinite extra dimen sion. 

Deffavet fc Menoul (|2007t ) discuss the possibil- 



distant GW sources appearing dimmer than they 
truly are, from the loss of GW energy flux to the 
extra-dimensional bulk. Indeed, in the presence 
of large distance leakage, flux conservation over a 
source-centered hypersphere requires that the GW 
amplitude scales with distance D from the source as 
h +x oc £)-C«K™-a)/2 j ( 3 ) 

where dim is the total number of space-time dimen- 
sions accessible to gravity modes. Thus, for dim > 
5, it deviates from the usual h +x (D) oc 1/D scal- 
ing. In principle, black hole merger events and asso- 
ciated host galaxies could thus reveal the leakage of 
gravity over scales of order a few Hubble distances, 
by comparison to purely electromagnetic Hubble di- 
agrams, which are immune to such leakage effects. 

This is only one of several possible modified 
gravity signatur es in the GWs from cosmologically- 
distant sources ([Deffavet fc Menoul 120071 ). Another 
class of signatures is related to the GW polarization 
signal, with possibly additional polarizations be- 
yond the two transverse q uadrupolar (+ X ) modes 
of general relativity (e.g., Will . 2006f h Signatures 
also exist in relation to the GW propagation ve- 
locity which, in modified gravity scenarios, can 
differ from the speed of light. In this respect, the 
possibility to time a cosmological GW, relative to 
an electromagnetic signal causally associated with 
the black hole merger, may offer unique diagnostics 
of large-scale modified gravity. This could reveal, 
for instance, that the phase of the GW signal de- 
viates from general relativistic expectations, once 
p ropagated over co smological distances. 

Kocsis et al. I (|2008l) have explored further the pos- 



ity that extra-dimensional leakage of gravity in 
DGP-like scenarios may lead to cosmologically- 



sibilities of measuring photon and graviton arrival 
times from a same cosmological source. A general 
difficulty with this approach is that there will be a 
systematic and a priori unknown delay in the emis- 
sion of photons, relative to the emission of gravi- 
tons, since the former must causally lag behind the 
perturbing gravitational event. This difficulty could 
be overcome if it were possible to calibrate the rel- 
ative timing of the photon and graviton signals at 
the source. 

Prior to coalescence, gas present in the near en- 
vironment of the black hole binary would be grav- 
itationally perturbed in such a way that it could 
radiate a variable electromagnetic signal with a 
period closely matching that of the leading-order 
quadrupolar perturbation induced by the coalesc- 
ing binary (see Fig. [5]). This would help identify 
the electromagnetic counterparts of specific GW 
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Fig. 5. [Left] From the numerical simulations of MacFadycn & Milosavlicvic (2008): snapshot of a gaseous disk gravitational- 
ly-perturbed by the time-dependent quadrupolar potential of a central binary. [Right] For such a system, in the inspiral phase, 
one may expect an electromagnetic (EM) source varying in brightness at a frequency approaching that of the GW signal re- 
lated to the quadrupolar perturbation. By matching the frequencies of the EM and GW signals, one could remove the delay 
in EM emission at the source, measure reliable offsets in arrival times between the two signals and possibly reveal drifts in 
frequency of the GW signal when it is propagated over cosmological distances. 



events. In addition, it may be possible to match the 
variability frequencies of the electromagnetic and 
GW signals. The offset in phase between the Fourier 
components of the two signals with similar frequen- 
cies could be used to effectively calibrate the intrin- 
sic delay in electromagnetic emission at the source. 
Late inspiral and coalescence can be tracked via the 
GW signal, so that the relative timing of the gravi- 
tational and electromagnetic signals may be known 
to within a fraction of the binary's orbital time. Any 
drift in arrival-time with frequency between the 
gravitational and electromagnetic chirping signals, 
as the source spans about a decade in GW frequency 
during the last 2 weeks before merger, could then be 
attributed to a fundamental difference in the way 
photons and gravitons propagate over cosmological 
distances. For instance, such a drift could occur if 
the graviton is massive, resulting in a frequency- 



symmetry could be violated in the gr avity sector 



e.g. 



Berti et al 



dependent propagation veloc i ty 
l2005t IWilll . l2006b iKocsis et all l2008h ^his tracking 
possibility is illustrated graphically in Fig. 

Interestingly, while Lorentz invariance has been 
extensively tested for standard model fields, Lorentz 



espec i ally on cosmological sc a les (e .g., ICsaki et al 



2001b IChung. Kolb & Riottol . l2002h . With a good 



enough understanding of the source, electromag- 
netic counterparts to black hole binary mergers may 
offer unique tests of Lorentz violations in the grav- 
ity sector, via the opportunity to match and track 
the gravitational and electromagnetic signals in fre- 
quency and phase. It may be possible, as the black 
hole binary decays toward final coalescence, span- 
ning a range of frequencies, to measure the delays 
in graviton vs. photon arrival times as a function 
of increasing frequency of the chirping signal. The 
consistency expected if Lorentz symmetry is satis- 
fied in the gravity sector could be tested explicitly 
for gravitons propagated over cosmological scales. 
To have any chance to perform such new tests of 
gravitational physics, one will need to identify the 
electromagnetic counterparts of coalescing pairs of 
massive black hole binaries as early as possible. This 
may be one of the strongest motivations behind 
ambitious efforts to localize these rare, transient 
events well before final coalescence. 
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4. New Possibilities with White Dwarf 
Inspirals 

As seen before, a large number of sources must be 
accumulated to turn a gravitational Hubble diagram 
into a high precision tool for cosmology. White dwarf 
frequently spiraling into (moderately) massive black 
holes may offer unique opportunities in this respect. 

An additional goal of the LISA mission is the 
detection of GWs emitted by compact objects 
being captured by massive black holes (the so- 
called Extreme Mass R atio Inspirals, or EMRIs). 
Barack fc Cutleil (2004) present a detailed parame- 
ter estimation analysis for this class of LISA events 
and show that they could be detected with good 
SNR out to a distance ~ 1 Gpc. In addition, these 
authors found that the sky localization errors for 
these events (~ 10~ 3 steradians at ~ 1 Gpc) are 
comparable to the case of black hole binary merger 
events. No electromagnetic counterpart is expected 
from EMRIs of dense neutron stars or stellar-mass 
black holes. On the other hand, partial (or total) 
disruptions of inspiraling low-density white dwarfs 
(WDs) could produce such counterparts and thus 
help identify host galaxies for such events. To the 
best of our knowledge, the possibility to find the 
electromagnetic counterparts of this subclass of 
LISA EMRIs and use them for precision cosmology 
has not been previously discussed in the literature. 

As a result of its relatively low density, a typical 
WD suffers complete disruption around a non- 
spinning black hole of mass < 10 5 M @ . This mass 
limit is increased to perha ps < 10 6 M Pl if the bla ck 



hole is spinning rapidly (|Ivanov fc Chernvakova , 
20061: iRathore et all . 120051 ). Even if the WD were 
not disrupted, some level of part ial shedding o f 



its lower density outer layers (e.g. Li et al. . 20021 ) 



may be expected as the inspiral proceeds. As the 
stream of debris from a partially stripped WD 
shocks against itself or the inspiraling WD (e.g., 
Rosswog et all [2008 ). some level of electromagnetic 
emission during such EMRIs is expected. A tidally- 
triggere d detonation of the WD is yet another pos- 
sibili ty ([Wilson fc Mathewsl . l2004t Dearborn et all 
2005). All these arguments point to the possibil- 
ity that electromagnetic counterparts may exist 
for a subset of all the EMRIs detected by LISA 
(the one corresponding to WD inspirals). The 
fraction of all EMRIs which involve WDs is, in- 



If electromagnetic counterparts for a subset of 
all LISA EMRIs can be detected, the rewards will 
be significant. Just like black hole binary merger 
events, EMRIs with uniquely identified counter- 
parts and host galaxies could be used to draw a 
gravitational Hubble diagram, with the significant 
advantages that, at redshifts as low as z < 0.2-0.3, 
weak-lensing due to line-of-sight inhomogeneities 
would be small or negligible (Fig. [2]) and that the 
D^-z relation at these low redshifts is strongly sen- 
sitive to the dominant dark energy content. Dark en- 
ergy becomes dynamically significant, affecting the 
expansion rate and geometry o f the universe, and 
modifying Dl, at z < 1 (see, e.g. lHuterer fc Turned . 



2001, for a review) . Indeed, as emphasized by 



Hu fc Haimal](|2003h . once cosmic microwave back- 
ground anisotropies are measured by Planck, Dl 
will be known accurately at z ~ 1000, and a low- 
redshift {z ~ 0) measurement will provide the best 
complement to constrain dark energy parameters. 

WD EMRIs detections with LISA will occur fre- 
quently, and they could perhaps be singled out on 
the basis of the comparatively low mass of the inspi- 
raling compact object. It remains to be determined 
how the LISA sky localization error evolves with pre- 
merger (or pre-disruption) time for such events and 
what is the nature of their electromagnetic counter- 
parts. But altogether, the various advantages that 
we have outlined point to the need for a detailed 
assessment of the potential use of WD EMRIs for 
precision cosmology with LISA. 



deed, expected to be substan tial (jGair et al.l . 12004 
Hopman fc Alexandeiir2006al lbh. 



5. Conclusion 

For centuries, astronomers have measured dis- 
tances exclusively with light. Direct gravitational 
measurements, gravitational Hubble diagrams and 
comparisons between the propagation of electro- 
magnetic and gravitational signals offer funda- 
mentally new ways to probe physics on cosmo- 
logical scales. The novelty involved in joint, time- 
constrained electromagnetic and gravitational mea- 
surements will require that special efforts be made 
to reach out across the GW and astronomy commu- 
nities. 

This work was supported by NASA grant 
NNX08AH35G. We thank A. MacFadyen for the 
permission to use one of his figures. 
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